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Chapter 1 

INTRODUCTION 


From engineering point of view any design should be safe and cost effec- 
tive. In civil engineering, during designing a system, most of the times 
attention is only paid to safety considerations resulting in uneconomic de- 
sign. On the other hand, to get a most economic i.e., optimum design 
many repeatative computations are to be carried out which are very much 
time consuming, tedious and monotonous. But the same can be achieved 
within a very short time with the help of a high speed digital computer. 
Sophisticated algorithms need to be used for the analysis to minimize the 
wastage of computer time. In the study reported here, the usefulness of 
such techniques in the optimal design of raft foundation has been pre- 
sented. 

For the economic design the objective function (cost function) is to be 
minimized subjected to some design restrictions arising out of material 
behavior and geometric considerations. As such, the problem is one of 
mathematical programming. Solution for such problems stretches the lim- 
its of conventional methods of calculus. Variational calculus is sometimes 
helpful but they are often limited in their range and versatility to tackle 
difficult problems. Many established optimization algorithms are available 
in literature. However, applications of these methods to new areas need 
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critical evaluation as no numerical method can be the panacea that can 
solve all problems and its efficiency is problem oriented. 

With the advent of high speed digital computers these methods started 
finding their applications in different branches of engineering. This led to 
the development of various nongradient and gradient based algorithms for 
unconstrained minimization. Some of these methods are quite stable and 
efficient in finding the minimum. Subsequently many direct and indirect 
methods of minimizing a function have been developed to solve constrained 
problems. With these developments more and more avenues were opened 
to tackle difficult systems which could not be taken up earlier. Appli- 
cations of these methods in Civil Engineering and specially in Structural 
Engineering began in the sixties. However, in Geotechnical Engineering 
their applications started only in the early seventies. There after for over 
two decades lot of these methods found their application in the area of 
slope stability analysis. However, their applications in the other areas of 
Geotechnical Engineering are very limited. So, there is a need to take 
up such studies. A brief review of the applications of these methods in 
geotechnical engineering is presented as follows. 

Design of a raft foundation involves proportioning and structural de- 
sign in which proportioning part mainly depends upon geotechnical as- 
pects, i.e., bearing capacity and settlement of the soil. Adequate bearing 
capacity and restricted settlement, these two are the primary and most 
important criteria for design of foundation. An intergated analysis of raft 
foundation from geotechnical and structural engineering aspects is very es- 
sential, which is rarely done. So there is a scope for optimum cost design 
of raft foundation considering these two aspects in design. 

Plenty of research work on the application of optimization methods 
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to be automated slope stability analysis are available in literature. Good 
reviews have been presented by various investigators (Bromhead, Groham, 
Nash, Choudhury, B hatt acharya) and, as such, those are not repeated 
here. Applications of these methods to other areas like optimal design of 
foundations (shallow or deep) and earth retaining structures are scanty. 
Some of these are as follows. 

Subbarao et al (1972) have developed a method for an optimal dimen- 
sioning of footings subjected to axial loading and uniaxial moment. But 
this method does not take care of the tilt of footing. Bhavikatti et al 
(1979) have proposed a method of optimal design of isolated column foot- 
ing using linear programming technique and reported about the progress 
of optimization and results of parametric studies. They concluded that 
using this method one can achieve 8 to 10 percent economy. But in this 
method settlement has not been considered. Madan Mohan ( 1984 ) has 
presented a methodology for settlement controlled optimum cost design of 
shallow footing. He worked on rectangular footing with generalized loading 
and used Sequential Unconstrained Minimization Technique. He has also 
reported a method of optimum plan dimensioning of a group of footings. 
Desai et al (1985) reported an optimum cost design of isolated footing and 
presented the results of different parametric studies. However, the study 
has considered only the axial loading on the footing. Mondal (1988) has 
studied the optimum cost design of shallow circular footings subjected to 
generalized loading. He has used Sequential Unconstrained Minimization 
Technique and reported about its success. Optimum design of a combined 
footing has been done by Pandian et al (1994). They have considered 
both bearing capacity criteria and settlement criteria to get the optimal 
solution for two columns. They restricted the settlement to 50 mm and 



4 


observed the changes of size of the footing with respect to the c han ge of 
column loading. Only width of the footing has been considered as a design 
variable. 

From the brief literature review presented above it can be concluded 
that till now no work has yet been done on the optimum design of raft 
foundation. For optimum cost design of such a foundation several trial 
designs need to be carried out satisfying all the structural and geotechnical 
design restriction. Even then enough cost reduction may not be achieved 
if the computations are done manually. All these can be easily achieved by 
resorting' to optimization formulation. As such, in this thesis optimal cost 
design of rigid raft with several column loads positioned in a symmetric 
way has been reported. 

Sequential Unconstrained Minimization Technique has been used in the 
presented study. Extended Penalty Function Technique (Kalvie and Moe, 
1971) supplemented by Powell’s Conjugate Direction algorithm (Powell, 
1964) coupled with quadratic interpolation technique (Fox, 1971) has been 
used. The above approach has been found to be quite efficient by various 
investigators in tackling problems even when the initial design vector is 
infeasible. As the success of any numerical scheme is problem oriented, 
it is necessary to have an in-depth study and a critical appraisal of the 
same to arrive at any meaningful conclusion regarding their efficacy and 
effectability. Such a study has been reported in this thesis. Parametric 
studies have also been carried out to bring out the effect of these these on 
the total cost. 

The organization of the thesis is as follows. 

In Chapter 1, a brief introduction and literature review, the motivation 
and scope of the work have been briefly outlined. 



Detailed analysis procedure identifying the objective function, the de- 
sign constraints, and nonlinear programming formulation to isolate the 
optimal solution is presented in Chapter 2. 

Results and discussions regarding the efficacy of the adopted technique 
and the influence of several design parameters on the overall cost have 
been presented in Chapter 3 - Finally in Chapter 3, the scope of future 
study has been reported. 



Chapter 2 


GENERALIZED PROCEDURE 
FOR OPTIMAL COST DESIGN 
OF RAFT FOUNDATION 

2.1 GENERAL 

Design of a rigid raft foundation involves, design of a slab, having a very- 
high thickness. In the conventional approach of designing, a rigid raft can 
be divided into a number of strips in both longitudinal and transverse 
directions. These strips are treated as combined footings subjected to 
several column loads. Structural design of the same is carried out using IS 
code of practice. 

Raft foundations can be adopted for steel as well as concrete super- 
structures. For concrete super-structure, columns can be supported di- 
rectly by the raft and no pedestal is required. But, for the steel structures 
pedestals are provided to connect steel super-structures with the concrete 
foundations. For maintenance purpose the joint should be placed at a 
higher level than ground surface. In this study for generalized column 
loading conditions both the options have been kept open for generalized 
column loading conditions. Limit State design method has been adopted 
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for concrete design. 

2.2 STATEMENT OF THE PROBLEM 

The Fig. 2.1 shows a rectangular raft with concrete columns while Fig. 

2.2 shows the same with steel columns. The connection of the raft with 
the super-structure may be either made of concrete or steel as shown in 
the figures. The problem is to find the optimal cost design of the raft for a 
given site condition, position of columns and the loading conditions (axial 
loading and biaxial moments as shown in Fig. 2. 3). The objective is to find 
the optimal cost design under the following two considered cases : 

1. Design of raft foundation for a set of concrete columns of given 
dimensions. 

2. Design of pedestals and raft foundation for a set of steel columns of 
given dimensions. 

2.3 ANALYSIS 

2.3.1 Assumptions 

In developing the generalized procedure of analysis the following assump- 
tions are made regarding the shape of the raft, columns, pedestals, spacing 
of the columns and the nature of the soil medium. 

1. Shape of the raft, columns and pedestals are rectangular or square. 

2. Spacings of the columns are same along x-direction. 

3. Spacings of the columns are same along y-direction. 


4. Soil is linearly elastic. 
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2.3.2 Design Variables 

When the raft has to be designed for concrete columns, number of design 
variables are 4. They are as follows : 

1. edge distance in x-direction, d x 

2. edge distance in y-direction, d 2 

3. thickness of the raft, d 3 

4. diameter of the bar, d 4 

When the raft has to be designed for steel columns, with the above 
variables three more variables, i.e., the dimensions of the pedestals 
are also to be included. They are as follows : 

5. width of the pedestal in x-direction, d 5 

6. width of the pedestal in y-direction, d e 

7. height of the pedestal, d 7 . 

2.3.3 Objective Function 

Total cost of the raft is considered as the objective function of the problem. 
It includes the cost of excavation, cost of construction, cost of concrete, 
cost of steel and cost of filling. If the raft is for steel super-structure cost 
of pedestal is to be included with the total cost of the raft. 

Objective function is given by, 

F — C exc + Cconc + C s t ee i + C fm 

— Rexc V exc 1” R-conc ^conc Rsteel 1'sieei Rfill ^ fill 
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where, C ex cj C conc ^ C steel and Cfni are the total costs of excavation, concrete 
(including construction charges), steel (including construction charges), 
and filling respectively. i? exc , R conc , R ste ei and Rjm are rates of the items 
and V exC i V conc , V s teei , Vfni are the volumes of excavation, concrete, steel and 
filling respectively. 

Total volume of excavation V exc is calculated in the following way, 

V exc f' R d br 


where, 

L and B are the length and width of the raft in m 
d br = depth of raft base from the ground level in m 

L = 2di + (n cp - 1) S x 
B = 2 d 2 + ( n rp - 1) S y 

where, di and d 2 are the design variables as described earlier 

n cp and n rp are the numbers of columns in each row parallel to x-direction 

and y-direction respectively. 

S s and S y are the spacing in m between two columns along x-direction and 
y-direction respectively 

Total volume of the concrete is calculated in the following way 

V CO nc R R ^3 V s teel 

where, 

d 3 = thickness of the raft in m 
other things are mentioned above. 

In case of steel super-structure with the above volume of concrete V conc , 
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the volume of concrete required for the pedestals is to be added. Concrete 
volume for pedestals is given by, 

Vped. ~ ft cp ftrp d§ dj 

As the volume of steel required for the pedestals is very less compared to 
the volume of concrete, volume of steel is neglected during the calculation 
of volume of concrete. 

Volume of steel is calculated as follows. At the points of maxi m u m 
bending moment occurring at the face of columns and in the span (the 
length between two columns) required reinforcement is calculated. The 
total volume of steel is given by, 

v ast = ±p sti 

1=1 

where, 

n = number of points where reinforcement is calculated 

Pst = reinforcement required at each point of calculation — N a l 

N = numbers of bar required 

a = area of cross-section of a single bar 

l = length of each bar taking the consideration of development length and 
anchorage length. 

As the thickness of a raft is very high, reinforcement is to be provided at 
both top and bottom faces. Method of reinforcement calculation is given 
later. 

Volume of filling V fiU is calculated in the following way. For concrete 
super-structure it is, 


V, 


fill 


Vexc kconc ^'ast ^'col 
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where, V co i —total volume of columns under the ground level up to top of 
the raft 

For steel super-structure it is, 

Vfill i'exc ^conc V ast ^'ped 

where, V ped = Total volume of pedestals under the ground level. 

2.3.4 Design Constraints 

For the design to be safe several design restrictions from the structural 
and geotechnical considerations are imposed. These are discussed under 
separate headings as folio wes. 

A. Structural Design Constraints 

1. The value of design variables d x and d 2 should be such that it 
will be at least equal to half of the column dimension in that direc- 
tion. Mathematically, 

di > w cx / 2 and d 2 > w cy /2 

in case of concrete super-structure. 

d\ > d 5 /2 and d 2 > d 6 /2 
in case of steel super-structure. 

where, and Wcy are column dimensions in x-direction and y-direction 
respectively. 

d 5 and d 6 are pedestal dimension in x-direction and y-direction respectively. 

2. If the dimensions of the design variables d x and d 2 are very large, 
then due to soil pressure high bending moment will develop at the outer 
face of col umn and it may cause failure of the structure. So the maximum 
dimensions of d x and d 2 are restricted to 1 metre. Mathematically 

di < 1 m and d 2 < 1 m 
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3. In case of rigid raft thickness of the raft is very high. But it should 
not be so thick requiring excessive amount of concrete and steel and pro- 
duce an uneconomic design. Again, it should not be so thin that the raft 
cannot be considered as a rigid one. So thickness of the raft is restricted 
as followes (Bowels, 1886) : 

d 3 > 0.75m and d 3 < 2.0m 

4. Diameter of the bar is also kept as a design variable. The size of the 
bar is restricted to : 

c 2 4 > 12mm and d 4 < d 3 /8 

In case of steel super-structure in addition to all the constraints already 
specified above the following constraints are applicable. 

5. Plan dimension of the pedestal at least 0.5m greater than the width 
of the column in the considered direction. Mathematically, 

do > Wcx + 0.5 and d 6 > Wcy + 0.5 

6. Maximum width of the pedestal is restricted to, 
d 5 < 1.5m and d 6 < 1.5 m 


7. Height of the pedestal is restricted to 
d 7 > 0.2m and 67 < 2.0m 

8. Height of the pedestal is restricted by the following manner also 
djr > dz + dr 


9. To check if the provided depth of the pedestal is sufficient or not, 
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for each pedestal the following relation should be satisfied [IS:456, clause 
33.1.3] : 


tana > 0.9 



+ 1 


where, 

q Q = calculated maximum bearing pressure at the base of the pedestal in 
N/ mm 2 

f C k = characteristic strength of concrete in N/mm 2 


Beside these structural constraints there are some other constraints 
applicable for detailing part of the reinforcement. As per IS:456 those 
constraints are given below. 

(a) Requirement of minimum steel [IS:456, Clause 25.5.1.1 (a)] : 

_ °' 85 
s^smin — r 

Jy 

where, f y = characteristic strength of steel in N/mm 2 

(b) Requirement of maximum steel [IS:456, Clause 25.5.1.1 (b)] : 

A S max ~ 0.04 b D 


where, 

b = width of the section 
D = overall depth of the section 

Provided steel must be in between the above limits, i.e., 


A sm in ^ A s t 5 ; A STnax 
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(c) Requirement of minimum spacing : 

^min ~ di 

where, S m i n = minimum spacing in between two bars 
d 4 = maximum diameter of bar used for design 


(d) Requirement of maximum spacing : 

Smax = 450 mm and S max = 3 d 3 
whichever is less. 

Provided spacing must be within the maximum and minimum values as 
specified above [IS:456, Clause 25.3.1], i.e., 


<S st 


< s„ 


. S st = spacing of the reinforcing steel. 

Reinforcement length should be such that it is in excess of the esti- 
mated length by a length called development length (L d ) estimated by the 
following relation [IS:456, Clause 25.2.1]. 


L d = 


d\(J s 

4 r bd 


where, d 4 = diameter of the bar 

a s = stress in bar at the section considered at design load 
r bd = design bond stress. 


For the members, subjected to combined axial load and biaxial bending, 
following test should be done [IS:456, Clause 38.6] : 



where, 

M ux and M uy are 


moments about x and y direction due to design loads. 
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M ux i and M ulgl are maximum uniaxial moment capacity for axial load P u , 
bending about x and y directions respectively. a n is related to 
Puz is given by : 

Puz = 0Aof c k A c + 0.75 fy .4 5C 

where, P u = axial load on the column 
A c = area of concrete 
A sc = area of steel. 

B. Geotechnical Constraints 

Constraints imposed on the design from geotechnical considerations are 

1. Settlement of the raft should be within permissible limits 

2. Imposed stresses should be less than the allowable bearing capacity. 
These constraints are described below. 

Settlements 

Main aspect of foundation design is to restrict the total and differen- 
tial settlements. Due to applied load on the raft, pressure on the soil will 
increase and it results some settlement of the underneath soil. Method 
used for calculation of stresses and settlement is as follows. For estimat- 
ing immediate settlement Boussinesq’s equation has been used. In case 
of cohesive soil there will be both immediate settlement and consolida- 
tion settlement and by adding up these two, total settlement is obtained. 
Consolidation settlement is estimated by using Skempton Bjerrum (1957) 
method. However, for frictional soils the settlement may be computed ei- 
ther by using the Boussinesq’s equation or by using S chmert mann method. 
The methods which are used are described below. 

Estimation of Stresses and Displacements 
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For calculating settlement estimation of stress increment at different 
points within the soil mass is necessary. As the raft is considered as a rigid 
one, stresses are calculated in the following way. Applied load on the raft 
may or may not be symmetric. So, with respect to the centroid of total 
area, eccentricity of the loading is to be determined. 


e x = 


E Pi Xi 
EP 


( 2 . 1 ) 


e 


y 


E Pi yi 
EP 


( 2 - 2 ) 


where, 

e x and e y are the eccentricities of the loading along x-direction and y- 
direction respectively. 

Pi = concentrated column loads 

Xi and yi are the distances of Pi from the centroid along x-direction and 
y-direction respectively. 


With these eccentricity at several points under the raft stresses are 
calculated using the following equation. 


EPi + R , M x + EPie y , M y + E,Pie x „ 
Pu = — _± T y± j i 


( 2 . 3 ) 


where, M x and M y are total concentrated column moments about x- 
direction and y-direction respectively. 

A = total area of the raft 

/ and I,, are moments of inertia of the raft about x-direction and y- 

x y 

direction respectively. 

x and y are the distances of the point from c.g. where stress is to be cal- 
culated along x-direction and y-direction respectively. 
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R = self weight of raft = L B W con 
W con = unit weight of concrete = 25 kN/m 3 

The whole raft is divided into several small rectangular elements. For 
each rectangular element stress is calculated. As the stress distribution 
under the raft is considered as linear, calculation of stress at four corners 
of each element is easy. From the stresses at four corners average stress 
at the centre of the element is calculated. If p a is the average stress at the 
centre of the element, 

Pa = Pa &£ 

is considered as the concentrated load at the centre of the element, where, 
a e is the area of a single element. At any depth z, stress increment is 
calculated using the Boussinesq’s equation given below. 


<*zi 


3&1 

2 7r r 5 


(2.4) 


where, 

cr zi = stress at z depth 

= point load applied at the ground surface 
Other quantities are as shown in the Fig. 2.4. 

At the four corners of the raft elastic settlement is calculated using Eq. 
2.5. Due a concentrated load P a acting at the surface of the soil settlement 
at any cylindrical radius r and at any depth z, the settlement is given by 
the following equation. 


pz - Pa 


1 + V 

2 TrE s r 


2(i + 


where, 

p z = vertical settlement of the soil 
E s = modulus of elasticity of the soil 


(2.5) 
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/tc. = poison’s ratio of the soil 

Here an option is kept to check the result using Mindlin’s approach. The 
equation is given below. Diagram is given in Fig. 2.4. 

n _ Pa 3 — 4/^ - | 8(1 — jU) 2 — (3 — , (z — c) 2 

16 7T G (1 [ R x ~R 2 + If 

(3 — 4 p)[z + c) 2 — 2cz 6cz(z + c) 2 ] 

M + M. . 

P a for each element will contribute to the settlement of each corner of 
the raft. 

For the ith element P ai causes p zi settlement then, total settlement at 
each corner can be stated as 


Pz — P zi 

With the same way settlement at the rest three corners can be calcu- 
lated. 

Consolidation Settlement 

If the raft has to be placed on fine grained soils, then there will be imme- 
diate settlement as well as consolidation settlement. Immediate or elastic 
settlement can be calculated either by using Boussinesq’s or Mindlin’s set- 
tlement equation (Poulos and Davis, 1974). Then consolidation settlement 
can be calculated using the following formula : 


S c = 


l + e c 


H log 


°o + Ap 


where, 

C c = compression index 

e 0 = initial void ratio of the soil 

H = thickness of the soil 


( 2 . 6 ) 



23 


o 0 = effective stress at the middle of the soil layer 

Ap stress increment due to applied load using Boussinesq equation 

If there are several layers, then total consolidation settlement can be 
calculated as : 

5 C = E^' 

where, 

S C j = settlement of the j th layer. 
j is the number of layers. 

Skempton-Bjerrum correction factor is applied on the result obtained 
from Terzaghi’s one-dimensional consolidation formula to take care of the 
three-dimensional consolidation effect. 


S actual — P S c 


where, ft = Skempton-Bjerrum correction factor 

The constant 0 is a function of dimension of the footing (breadth, B), 
depth of the compressibe layer (z) and porepressure coefficient ‘A’. The 
variation of 0 with ‘A’ for various values of zj B and for circular/strip 
footings are produced in graphical form. 

As per IS:1904-1978, maximum allowable settlement for the raft placed 
on either clay or sand is 100mm and 50mm respectively. Differential settle- 
ment is calculated by estimating the difference of total settlement at two 
extreme corners. If pi is the settlement at corner 1 and P 2 Is the settle- 
ment at corner two and l is the distance between those points then tilt (0) 
is given by 


l 
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Maximum tilt is restricted to 0.003, i.e., 

B < 0.003 


Schmertmann Method 

If CPT data is available for the site then this method can be used 
to determine total settlement. In this method the strain is assumed to 
increase from a minimum at the base of to a maximum of 0.6 at S/2 depth, 
then decrease and reaches zero at the depth equal to 2 B. The area of the 
strain triangle gives settlement. The triangular vertical strain influence 
factor distribution diagram as proposed by Schmertmann (1978) is shown 
in Fig. 2.5. 

The equation for calculating settlement is : 

2 B j 

Az (2.7) 

0 

where, 

S = total settlement 

q n = net foundation base pressure = (q - q 0 ) 
q = total foundation pressure 

q' o = effective overburden pressure at foundation level 
A z — thickness of the elemental layer 
I z = vertical strain influence factor 
Ci = depth correction factor 
C 2 = creep factor 

The equations for Ci and C 2 are : 

I 

Ci = 1 - 0.5 (^) 

Qn 
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C 2 - 1+ 0.2 l O g l0 (~) 

where, 

t =time in years. 

E s can be calculated from the following equation : 

E s = aq c ^ ( 2 . 8 ) 

a depends upon the type of the soil. 
q c = cone resistance A computer program has been developed to estimate 
the factor I z corresponding to any depth (< 25). 

Bearing Capacity 

This is the other important criteria for the design of a foundation. If 
the soil pressure exceed allowable bearing capacity of the soil then there 
will be shear failure of soil underneath the raft. Bearing capacity of the 
soil is calculated by using Hansen’s ultimate bearing capacity equations 
(Bowels, 1982, 1988). 

For c — 4> — soil : 

Quit = c N c s c d c i c + q N q s q d q i q + ^B 7 iV 7 s 7 d 7 i r (2.9) 

For <fi — soil : 

quit — q N q s g d q i q + - B 'y N 7 s 7 d 7 i 7 (2-10) 

For c — soil only : 

q u it — 5.14 Su (1 + s c + d c )+q (2-11) 

where, bearing capacity factors 


N q = e ir un ^tan 2 (45-f <£/2) 



N c = (N q — 1) COt 0 


^7 — l-5(iV ? — 1 ) tan<£ 


shape factors 

s' c = 0.2 B/L 


<s c — 1 + Ng B/L 
s q = 1 + (B/L) taiuj) 
s-y = 1 - 0.4 B/L 


depth factors 


dq 


for all (j> 



0.4 D/B if D<B 

0.4 tan ’ 1 D/B if D > B 


d>c — 


1 + 0.4 D/B if D < £f 

1 + 0.4 tan” 1 D/B if D > B 


f 1 + 2 tan cj)(l - sin ) 2 D/B if D < B 

| 1 + 2 tan <f>(l - sin ) 2 tan -1 D/B if D < B 


dy = 1.0 


As all the loading conditions are considered as vertical, 
factors are taken as 1, i.e., 


%e — Iq lny 1.0 

. S u = undrained shear strength 
c = cohesion intercept 
<j i ) = angle of shearing resistance 
D = depth of the foundation 
B = width of the foundation 
9 = 7 D 
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Inclination 
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In Hansen s equations other factors like, ground factor and base factors 
are also there. But here, ground and base of the raft are considered as 
horizontal and as such those factors are not considered. 

Factor of safety for raft against bearing capacity failure varies from 1.7 
to 2.5 (Bowels, 1982, 1988). 

So, allowable bearing capacity is calculated by the following way, 

Qall = Qult/F.S. 

Here the adopted value of factor of safety is 2.0. 

In ultimate bearing capacity equations for and c - <p soil there is a 
term including the dimension of the raft (width of the raft B). It seems 
that if B increases indefinitely then ultimate bearing capacity will also 
increase indefinitely to an infinite extent. But actually it does not occur 
in that way. Up to a certain size of foundation q u i t increases and then it 
becomes constant though size of the foundation increases. Bowels (Bowels, 
1988) has suggested a reduction factor r 7 which can be calculated using 
the following formula if B > 2m : 

r 7 = 1 — 0.25 log(B/K ) 


K = 2.0 for SI units. 

This reduction factor is used in the following way, 

0.5 7 B iV 7 s 7 d 7 r 7 

This takes care of the effect of the size of the raft beyond a certain limit. 
Now the bearing capacity constraints can be stated as : 

If q r is the produced stress at the four corner of the raft then q r < q a u 


Check for Relative Rigidity of raft 
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It is well known that for design of beams on elastic foundation if the 
maximum spacing between two columns (S) is such so that the product of 
A and S is less than 7r/4 then the beam is called as rigid onej A is a factor 
having the dimension of 1/ length and can be calculated using the following 
formula : 


A = 


K t B s 
4 E b I b 


( 2 . 12 ) 


where, 

K s = modulus of subgrade reaction of the soil 
B s = width of the strip 

E b = elasticity modulus of the material of the strip (concrete) 
= 5700 SU 


I b = moment of inertia of the section of the strip = 1/12 B d z 
d = depth of the beam 

Modulus of subgrade reaction K s can be calculated using the following 
equation (Vesic, 1961a, 1961b) : 


K s = 0.65 



E, 

l-/i 2 


( 2 . 13 ) 


where, 

E s = elasticity modulus of the soil 
fj, = poisons ratio of the soil 


For all the strips dimension of the width should be less such that the 
value of A should be sufficiently less so that the product of maximum 
spacing of the column in the raft and A will be less than 7r/4. In that 
case only the adopted method is applicable and as such, this check is very 


essential. 
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2.4 DESIGN 


Once the different dimensions of the raft, depth of the foundation, load- 
ings and soil parameters are known the structural design can be done. The 
method adopted here to design a raft is the conventional strip method. In 
this method the whole raft is taken as a combination of several strips ac- 
cording to the number of column rows in both x and y directions. Then 
each strip is designed separately. For design the following steps are fol- 
lowed. 

1. All the vertical column loads and concentrated column moments 
coming on the whole raft about y-direction is taken into account. The 
maximum and minimum pressure at the two ends of the raft along x- 
direction is estimated by using the following equation. 


PX 1,2 = 


X) P, + .Rs 

as 


My + Pfix 
— L _7 ~ 2/ 


( 2 . 14 ) 


where R s = total weight of the strip 
i = number of columns on that strip 

2. The raft is divided into several strips, according to the number of 
column rows, parallel to the x-direction. 

3. A single strip is considered parallel to x-direction. Eccentricity of 
the loading is calculated using Eq. 2.1. For this eccentricity and the 
concentrated moment on that strip maximum and minimum pressure is 


calculated at the two ends of the strip using Eq. 2.14. 
Where, i = number of columns on that particular strip 


Rx= R 3 = self weight of the strip. 

The stresses are psi, 2 - I* 1 Fig* 2.6 calculation of p^i ,2 and psi ,2 are shown. 

4. The average of pr,, and ps,, 2 is calculated and the average stress is 


pa i, 2 * 
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5. From the average stresses pa 1)2 total soil pressure S P is calculated by 
the following way : 

Sp = 1/2 (pa 1 +pa 2 ) L s B s 

where, L s = length of the strip 
B s = width of the strip 



(b) 


Figure 2.6: Stress calculation for design, (a) Considering whole raft, (b) Considering 
one strip 


6. If the total load on the strip is L T including the self weight of the 
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raft then the load reduction factor Lrf can be calculated in the following 
way. 

Lrf = Sp/Lr 

• All the concentrated column loads and the self weight of the strip is 
multiplied by L RF and the obtained loads are used to calculate design 
moments. 

This has been done to satisfy force equilibrium condition for that par- 
ticular strip (Bowels, 1982). The same thing has been done for moments 
for that particular strip. If the moment reduction factor is given by Mrf-, 
then all the concentrated moments are multiplied by Mrf and the ob- 
tained moments are used subsequently in the design moment calculation. 
At each face of the column and at the span in between two columns bend- 
ing moments are calculated. At the face of the column maximum bending 
moment is considered as design bending moment at that column point and 
maximum span moment is considered as design span moment at that span. 
The following equation given in Appendix E of IS:456 is used. 

the required reinforcement at those points. 

M u = 0.87 f y A st d (1 - yLjt) ( 2 . 15 ) 

where, 

M u = moment for which reinforcement has to be provided 

f y = characteristic strength of steel in Nj mm 2 

f ck = characteristic strength of concrete in N/mm 2 

d — effective depth of the section 

b = width of the section 

A st — area of reinforcement 

If the bending moment is sagging at the column points reinforcement is 
to be provided at bottom and if the bending moment is hogging at the span 
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reinforcements are to be provided at the top. But in both cases minimum 
reinforcement is to be provided at the opposite faces as the thickness of 
the raft is very high. 

With the same procedure required reinforcement can be calculated for 
all other strips parallel to x-direction. To calculate the reinforcement for 
the strips parallel to y-direction all vertical column loads and the concen- 
trated column moments about x-direction are to be considered. Other 
things for calculation are same as given in the steps stated above. 

In case of steel super-structure during design of raft the column loads 
are increased by ten percent to take care of the weight of the pedestal. 
Other design procedure is same as stated above. 

Check Against Shear 

For each column which \§ supported by the raft two way shear, i.e., 
punching shear test is must. Critical section for punching shear is at a 
distance of d 3 /2 from the face of the column. Here is the effective depth 
of the raft. Shear stress can be calculated by 

Pi 

Teal — , • 

lpe “3 

where, 

Pi = vertical load on i t h column 

l pe = periphery of the area at a distance of (d 3 /2) from the column face. 
For interior columns 

lpe = 2 (wex + Wcy + 2 d 3 ) (2.16) 

For side columns, if d 3 / 2 is greater than the edge distance from the face of 
the column then 

l pe = 2(w cx + d 3 ) + ( w C y 4 - d z /2) 
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or 

Ipe = ~{ w cy + <4) + (U-'cx T ($ 3 / 2 ) 

Otherwise, Eq. 2.16 is to be used. For corner columns, if <4/2 is greater 
than the edge distance from the two perpendicular faces of a corner of a 
column then 

Ipe = (Wcx + IL'cy + 4 ) 


Otherwise, Eq. 2.16 is to be used. This must not exceed k s r c [IS:456, 
Clause 30.6.3.1]. where, 

k s = (0.5 + fi c ) 
k s 1.0 

p _ short side o f the column 
long side of the column 

r c = 0.25v7cit for Limit State method of design. All columns have to satisfy 
this condition. 

2 . 5 FORMULATION OF MATHEMATICAL PRO- 
GRAMMING 

Determination of minimum value of objective function subject to con- 
straints as stated earlier can be stated as a generalized optimization prob- 
lem as follows : 

Find D m such that 
F(D m ) is minimum 
subject to <?j(D m ) < 0.0 

where, D m is a function of design variables 
< 7 j(D m ) are constraints. 
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The problem is converted into an unconstrained optimization problem 
with the help of Extended Penalty Function Method (Kavlie and Moe, 
1971). The Sequential Unconstrained Minimization of the developed com- 
posite function is then carried out in conjunction with Powell’s Conjugate 
Direction Method for multidimensional search with Quadratic Interpola- 
tion Technique for linear minimization. The composite function (f)(D, r*) 
developed by blending the objective function and constraints is as follows 

^(D,r i ) = F(D)+r i £G[ a (D)) 

j=l 

The function G[<?j(D)] is chosen as 


G[ gj m = 


where, e = —rk/S 


1/&(D) when g s { D) < e 
2 - ~{ - - D when ffj(D) > e 


Tk = penalty parameter and <5 = transition parameter. 

Advantage of this method is that even if the initial design vector is in the 
infeasible region the algorithm is capable of finding the feasible optimum 
solution. For various engineering problems most of the times initial feasible 
design vector is difficult to guess. In such cases, the developed technique 
is expected to be helpful in finding a feasible solution. 


2.5.1 Convergence Criteria 

If the relative change of objective function and design variables between 
two consecutive cycles satisfy the following conditions then the convergence 
is assumed to be achieved. 


Fk — Fk- i 

F k 


< e t 


. (2.17) 
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k = number of iterations 
F = objective function. 
e, = very small numerical quantity. 



Chapter 3 

RESULTS AND DISCUSSIONS 


A program has been developed in FORTRAN language to obtain the opti- 
mal solution following the procedure outlined in the previous chapter- All 
computations have been carried out using Pentium processor. Effective- 
ness of the suggested method has been demonstrated through parametric 
study using different soil parameters (like c, o, E s » V, C c /(1 + e a ) ). 

To estimate the total cost of the raft the rates used for concrete, steel 
and other construction items are according to the Central Public Works 
Department’s schedule rate for the year 1997. The rates are presented in 
the Table 3.1. The characteristic strength of concrete and steel (f c k and 
f y ) are taken as 20 N/mm 2 and 415 N/mm 2 respectively. 

Using the developed program it is possible to design a raft on any type 
of soil. The soil may be purely cohesive, purely frictional or a general 
cohesive-frictional soil. It can also be used when the soil data is available 
in terms of CPT value with depth for the considered site. The permis- 
sibe value of settlement has been taken to be 50mm if it is not specified 
otherwise. 

For a given set of loading shown in Figure 3.1 parametric study has 
been carried out and the obtained results are as follows. 

The developed computer program has been varified by computing the 
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Table 3.1: Rate of the items in Rs/m 3 


Particulars | rates j 

cost of excavation 

15 

cost of concrete 

1230 

cost of steel 

1670 

cost of filling 

8 


obtained solution for different strips with manually computed values. The 
problem and the results obtained by computer program and by the hand 
calculation are shown in the Table 3.2. The results show good agreement 
with each other. 

The influence of the initial design vector on the optimal solution has 
been studied. The results are presented in Table 3.3. Soil condition is 
same as given in Table 3.2. It can be seen from the table that starting 
from widely differing initial design vectors the optimal vectors is practi- 
cally identical and the costs differ only marginally (0.2 %). So, the optimal 
solution can be considered to be the global solution. Convergence of the 
numerical scheme with decreasing sequence of the penalty parameter r k is 
shown in Fig. 3.2. It can be seen that as r k decreases the composite func- 
tion and the objective function tend to converge for r, less or equal to 20 
showing that minimization of the composite function led to the minimiza- 
tion of the objective function. From Fig. 3.3 it is seen that as the number 
of function evaluations exceeds about 350 both the composite function and 
objective function converge. It is observed from Fig. 3.4 that when the 
number of iterative cycles exceeds 2 the convergence is gradual and ex- 




39 


Table 3.2: Comparison between hand-calculation and program-obtained value 



Table 3.3: Effect of initial variables on the result 


1 Initial variables 

(m) 

Final variables (m) 

Cost (Rs) 

<*1 

d 2 ‘b 

d 4 

dl 

d 2 

b 

d 4 


1 

1 1 

0.02 

0.15 

0.15 

0.75 

0.012 

48577 

4 

4 2 

0.6 

0.154 

0.150 

0.75 

0.0132 

48680 





Column no. 

P 

Mx 

My 

1 

250 

-40 

40 

2 

250 

-20 

30 

3 

300 

-40 

30 

4 

200 

-50 

30 

5 

200 

-20 

40 

6 

250 

-40 

30 

7 

200 

-30 

40 

8 

250 

-40 

20 

9 

250 

-20 

40 


Figure 3.1: Sample loading from the super-structure 
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cellent. These figures show that the adopted method is quite effective in 
isolating the optimal solution. 



Figure 3.2: Convergence of objective function and composite function with respect to 
Penalty parameter 

The stress and settlement calculations have been carried out by dis- 

V 

cretizing the plan area of the raft into several small elements and taking 
the effect of each element at each corner of the raft. Equations 2.4 and 
2.5 have been used to calculate stress and settlement. As the application 
of Boussinesq’s equation cannot take into account the effect of the embed- 
ment of the foundation in the soil medium, the results have been checked 
by using Mindlin’s solution. The comparative results are shown in Table 
3.4. It can be seen from the table that the results do not vary much (0.05 
% only). For all the cases value of E s is 50000fciV/m 2 . This shows that for 
the given problem the embedment of the foundation does not play any 
significant role on the overall cost and so also the other parameters like /i, 
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Figure 3.3: Convergence of objective function and composite function with respect to 
Number of function evaluations 



Figure 3.4: Convergence of objective function and composite function with respect to 
number iterative cycles 





43 


Table 3.4: Effect of Boussinesq’s and Mindlin’s solution, on the results 


Smax 
in mm 

c(kN/m 2 ) 

0 

(in degree) 



cost in Rs 
by Boussinesq 

cost in Rs 
by Mindlin 

ElHH 


10 

0.45 

0.06 

48557 

48556 1 

|TT3iyp 


15 

0.45 

0.06 

48545 

48586 




0.45 

0.06 

48524 

■EH91 


10 

30 


0.05 

48557 

MiliM 




0.3 

0.05 

48521 

48522 



| 

0.3 

0.05 

flBENI 

48565 

mm 




0.05 


48565 


C c /{ 1 + e o)i <t>, c and settlement limit. This is due to the fact that for these 
paramaters settlement constraints are inactive. 

A study is made to estimate the ammount of saving that can be achieved 
over a manual design. For this a few trials are made and the best out 
of these is fed to the computer as the initial design vector. At the first 
iteration it gives similar results to the values obtained by hand calculation. 
After 9 iterations final optimal results are obtained. The results are shown 
in Table 3.5. It can be seen that at the optimal point the settlement 
constraints are all active and so also some geometric and bearing capacity 
constraints. A net saving of Rsl 1,892/- can be obtained which is about 
19.65 % of the best possible cost that has been acieved manually. This 
demonstrates the necessity of adopting optimization procedure to get the 
mimimum cost design. 

Parametric Study 

For saturated fine grained soils under undrained condition 0=0), the 
effect of the undrained shear strength ( S u ) on the overall cost of the raft has 
been studied by varying the S u value ranging from 20 kPa to 40 kPa. The 
results are presented in Fig. 3.5. It is observed from the figure that when 
























Table 3.o: Comparison between best hand-calculated result and optimum result for 
same soil and loading condition 


Soil condition : 

_ = 0.499 Cc/d+Cy) =0.07 

L = lOkt’o GWT at 0.0 m depth 

Y =20kN/m 3 

Hr = 25 MPa 


Initial variables ( d j,d ^ d 3* d 4) 

0.5, 0.5, 0.75, 0.016 
Initial area ratio (Ar) 

25.724 

Initial constraints ( first 8 constraints are geometric constraints, next 4 arc differential 
settlement constraints, next 8 are bearing capacity constraints and last 4 are total 
settlement constraints) 

-0.5, -0.5, 0.0, -1.25, -0.35,-0.35, -0.0040, -0.0777 
0.0034 38, -0.0030, -0.00331, -0.0033, -659. .933, 699,816, 

-645.239, -645.239, -134.26, -94.3768, -188.837, -148.954, -0.014, 

-0.0155, -0.0125, -0.0109 
Initial cost 
Rs 60518 


Final variables ( d j. d?*d 3, d 4) 

0.151, 0.155, 0.75, 0.012 
Final area ratio(Ar) 

1.806 

Final constraints ( first 8 constraints arc geometric constraints, next 4 arc differential 
settlement constraints, next 8 are bearing capacity constraints and last 4 are total 
settlement constraints) 

-0.8482, -0.8447, -0.000054, -1.2499, -0.00174, -0.005277, 

-0.000411, -0.0813, -0.00354. -0.00294, -0.00334, -0.003421, 

-722.765, -777.384, -648.274, -702.893, -160.154, -105.535, 

-234.645, -180.025, -0.00965, -0.00115, -0.00774, -0.00655 
Final cost 
Rs 48626 
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S u > 30 kPa the undrained shear strength value does not have any influence 
on the cost of the raft. The influence of the undrained shear strength value 
on the cost is appreciable for strength value less than 30 kPa. As the the 
strength value is decreased from 30 kPa to 20 kPa , the rate of increase 
in the cost increases at an increasing rate. Rate of increase in cost is 9 
% when S u is decreased from 30 kPa to 25 kPa (i.e., 16.7 % decrease in 
strength ) and 20 % when it is further decreased to 20 kPa from 25 kPa (20 
% decrease in strength). This increase in the cost when the strength value 
decreases can be justified from the study presented in the Fig. 3.6 where 
the influence of the varying strength value on the area of the raft has been 
shown. To do this, a parameter ‘area ratio’ (A r ) has been introduced 

, t 1 

64000.- <+> = 0 

Es = 40 MPa 

62000 - \ u = 0.45 

Cc/(l+e 0 ) = 0.05 

60000 - •. GWT at 0.0 m depth " 

„ 58000 - \ ?= 20 kN/ ” 3 

•9 56000 - 

4-J 

C/3 

O 54000 - 

+.. 

52000 - ■•••. 

50000 - 

*+ 

48000 t I i I I 

20 25 30 35 40 

S u in kPa 

Figure 3.5: Variation of cost with S u 
which is given by the following formula: 



+. 


_j _ 


<{> = 0 

Es = 40 MPa 
u = 0.45 

Cc/(l+e 0 ) = 0.05 
GWT at 0.0 m depth 
Y= 20 kN/m 3 
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where, 

A a = base area of the super structure 
.4x = area of the raft 

It can be seen from the figure that as strength value decreases the change 
in A r is not significant till the strength value does not fall below 30 kPa; 
below this value the increase in the ‘area ratio’ is very sharp for decreasing 
values of S u . It can be seen that the increase in ‘area ratio’ is 9 % when 
strength value decreases from 30 kPa to 25 kPa while the same is 22 % 
when the strength value is decreased to 20 kPa from 25 kPa. As S u is 
decreased below 30 kPa the bearing capacity value decreases resulting in 
the increase area of the raft and consequent increase in the cost- When S u 
is greater than 30 kPa , the bearing capacity does not control the design. It 
is interesting to note that over the whole range of S u values varying from 
20 kPa to 40 kPa the depth of the raft remains almost constant throughout 
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(Fig. 3.7). This signifies that the increase in cost is only due to the 
increase in plan area of the raft. From Fig. 3.8 one can observe that the 
cost of concrete predominantly influence the cost of the raft. Over a wide 
range of undrained shear strength the ratio of the cost of concrete and 
cost of steel almost remain same and the ratio is about 5.2 to 5.3. So to 
save the cost one can design the raft as a flexible one which will give a 
low value thickness resulting in ultimate savings in the cost of concrete. 
However, in such a case the effect of differential settlement at the different 
column points have to be taken into account on the overall behaviour of 
the structure. Therefore, one must adopt integrated design procedure 
considering the super-structure and foundation as a single uni t. 



Figure 3.7: Variation of thickness of the raft (d 3 ) with respect to S u 


From Fig. 3.9 it can be seen that for purely cohesionless soils for given 
values of </> and other parameters, the variation of p over a range of 0.35 to 
0.499 does not play any significant role on the cost as well as on the optimal 
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dimensions as the design is not controlled by settlement. Settlement is 
always less than the permissible value. But, from Fig. 3.10 it can be 
noticed that for E s below 25 MPa effect of E$ on the raft is greater. At 
this stage design is strictly controlled by settlement and plan area of the 
raft is highly increased resulting a high ‘area ratio’ ( A r ). But, the depth of 
the raft did not change in any case. This is due to the fact that thickness 
required from shear and bending moment criteria is much less than the 
geometric constraint (dz > 0.75m) and, as such, this constraint decomes 
active during optimization. For the sake of brevity all the results are not 
presented here. 

From Fig. 3.11 it can be seen that for the given set of soil parameters, 
C c /( 1 + e 0 ) does not have any significant effect on the cost, till it reaches 
the limiting value of settlement up to which the design is not governed by 
the settlement criteria. But, as the value of CJ{ 1 + e 0 ) increases the rate 




Cost in Rs Cost in Rs 
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50000 

49500 

49000 


48500 

48000 

0.34 0.36 0.38 0.4 0.42 0.44 0.46 0.48 0.5 

Poisson’s ratio , jx 

Figure 3.9: Variation of cost with respect to fi for different values of E s 

64000 
62000 
60000 
58000 
56000 
54000 
52000 
50000 
48000 

Figure 



20 25 30 35 40 45 50 


Es in MPa 

3.10: Variation of cost with respect to E s for different values of <f> 
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of settlement also increases and when it reaches the limiting value ( 50 mm) 
the design is fully controlled by settlement. The variation of obtained ‘area 
ratio’ due to the change in C c /{ l+e 0 ) are shown in Fig. 3.12. It is seen from 
this figure that for E s equal to 25 MPa the result is not affected much and 
the settlement reaches nearer to the limiting value (50mm) as the design is 
then fully controlled by settlement. In Fig. 3.12, for E s equal to 25 MPa 
the result is not affected much till the value Cc/( 1 + e 0 ) reaches 0.105 but 
when it is 0.114 the cost of the raft is very high. Similarly, for E s equal 
to 50 MPa , up to the value of C c /{ 1 + e 0 ) equal to 0.15 the cost of the raft 
does not change. When it is increased to 0.154 cost is increased by about 
three thousand rupees. But, the cost is maximum when it reaches 0.164 
beyond which no solution could be formed. But, there is no change in the 
obtained thickness of the raft, which is only 0.75m. 



C c /(l+e6) values 

Figure 3.11: Variation of cost with respect to C c /{ 1 + e 0 ) for different values of E, 
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C c /(l+eo) values 

Figure 3.12: Variation of ‘area ratio’ (A r ) with respect to C c /( 1 + e 0 ) for different 
values of E s 

To see the effect of permissible-limit of settlement on the optimum cost, 
four limiting values of settlement i.e., 25mm, 50mm, 75 mm and 100mm have 
been chosen. For achieving enhaunced settlement the column loads are 
increased sequentially by some percentage. The results are presented in 
Table 3 . 6 . It is seen that only for 5 % increase in the initial loading as given 
in Fig. 3 . 1 , there is a great increase in the cost when maximum settlement 
is limited to 25 mm. At the same loading condition if the settlement limit 
is increased to 50 mm and 75 mm the optimal cost is reduced to a minimum 
level. No further change in permissible limit of settlement causes any 
change in the cost. Similar behaviour is observed when the permissible 
settlement limit is increased to 75 mm from 50mm and the raft can take 280 
% of the initial loading. When the settlement limit is increased to 100mm 
from 75 mm, the raft can carry 450 % of the initial loading. Table 3.6 shows 




the effect of maximum settlement limit on the optimum cost as well as on 


1 able 3.0: Results for different maximum limiting settlements 


Soil Condition 
c = lOki’a 
<j> = 15° 
p =0.499 

Es = 50 MPa y = 20 kN/m 3 

C c /( 1 +e $ = 0.07 GWt at 0.0 m depth 

Optimum Results for Maximum Settlement Limits 

Loading (muiti- 
-ple of loading 
given in 

Fig. 3.1) 

Smax = 

25 mm 

Smax = 

50 mm 

Smax = 

75 mm 

Smax = 

100 mm 

1.05 

cost Rs 50937 

Ar = 6.53 
d 3 = 0.75 m 
d^= 12 mm 

cost Rs 48591 

Ar= 1.576 
d 3 = 0.75 m 
d 4 = 12 mm 

cost Rs 48555 

Cost Rs 48555 

2.8 

no solution 

cost Rs 7 1 027 

Ar = 46.903 
d 3 = 0.75 m 

d 4 = 12 mm 

cost Rs 48627 

Ar = 1.592 
d 3 = 0.75 m 
d 4 = 12 mm 

cost Rs 48661 

4.5 

no solution 

no solution 

| 

cost Rs 73210 

Ar = 53.50 
d 3 = 0.75 m 
d 4 = 12 mm 

cost Rs 48815 

Ar = 1 542 
d 3 = 0.75 m 
d^= 12 mm 

6.5 

no solution 

no solution 

no solution 

cost Rs 74130 

Ar = 55.49 
d 3 = 0.75 m 

d4= 14 mm 


the ‘area ratio’ (.4 r ). Here also, the thickness of the foundation does not 
affect the cost and the value always being at its lower limit, i.e., 0.75m. In 
all the cases, the high cost, can be attributed to the setlement limit and 

corresponding increased ‘area ratio’ (.4 r ). 

program can also he used if the raft is placed soil strata for which 
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CPT data is available. To demonstrate its use, the test data reported by 
Singh (1982) for IIT Kanpur Campus has been used (Fig. 3.13). The value 



Figure 3.13: Variation of cone resistance with depth 

of a for Campus soil is 9 (Prakash Rao, 1985). Using Equation 2.8, E s is 
calculated for different q c and using Schmertmann Method total settlement 
is calculated. Obtained optimal solution is given in Table 3.7. 

Table 3.7: Results for layered cohesionless soil 


optimum dimensions 
in m 

optimum cost 
in Rupees 

area ratio 

6.3 X 6.3 X 0.75 

Final vector 

0.15, 0.15, 0.75, 0.012 

48526 

1.581 


In Fig. 3.14 reinforcement detailing for the above case is shown. 
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To see the effect of loading on the dimension of the raft, the load on 
the central column is increased and the obtained value of ‘area ratio’ A r is 



0 500 1000 1500 2000 2500 3000 


Load at central column.in kN 

Figure 3.15: Effect of increased loading on ‘area ratio’ A r , for c-soil 

can be observed from the figures that for (j> = 0 case the effect of increase 
in load on the central column is more pronounced than that in c — cf> -soils. 
The thickness of the raft is again the same ( 0.75m ) for all the cases and 
equal to the lower limit imposed on the depth of the raft. 

For the same set of loading the effect of spacing of the columns on the 
area of the raft has also been studied. It can be seen from Fig. 3.17 that for 
<j) = 0 case when the spacing is less than 3m the effect is highly pronounced 
and the ‘area ratio’ (A r ) suddenly increases. This is due to the bearing 
capacity requirement. If the spacing is less then total base area of the 
super-structure decreases, and to get required bearing capacity area of the 
raft increases. But, for general c — <j> -soil A r does not vary much with the 
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Figure 3.16: Effect of increased loading on ‘area ratio’ A r , for c- <p -soil 


increase of the spacing. It is shown in Fig. 3.18. However this observation 
may be problem specific. 

The results discussed till now are for designing the raft for concrete 
super-structure. In the program there is a provision to design a raft for 
the steel super-structure. For the same loading and for only cohesive soil 
a sample problem has been solved. In Table 3.7 comparison between the 
results obtained for steel super-structure and concrete super-structure has 
been produced for same soil condition as shown. It can be observed from 
the results that cost for steel structure is more and the additional cost can 
be attributed to the cost of the pedestals. 




5 ' 
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Table 3.8: Comparative results for steel and concrete super-structure 


Soil conditions : 

<t> = 0 GWT at 0.0 m depth 

c = 20 kPa 

Es = 40 MPa 

H = 0.45 

Cc/(l+e ) = 0.05 

y = 20 kN/m 
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3.1 CONCLUSIONS 

On the basis of obtained results and discussion presented earlier following 
conclusions can be drawn. 

1. The developed procedure where in the optimal cost design of a 
rigid raft has been formulated as a nonlinear programming problem has 
been found to be quite efficient in isolating the optimal solution. The 
composite function developed by using extended penalty function approach 
proposed by Kavlie and Moe (1971) behaved very nicely during sequential 
minimization and converged to the optimal solution solution very fast. 
Adoption of Powell’s Conjugate Direction method for multi-dimensional 
search and quadratic fit for linear search has been appropriate as evident 
from the fact that during the search for the minimum the minimization of 
the composite function led to the minimization of the objective function 
without any difficulty. 

2. For saturated fine grained soils under undrained condition, the 
undrained shear strength plays a great role on the optimal dimensions 
of the raft when its value is less than a limiting value. However, above 
that value the effect of undrained shear strength is negligible. Thickness 
of the raft is not affected by the variation in the undrained shear strength. 

3. For purely frictional soils there is no effect of <f> (20° <(f>< 30°) on the 
optimum cost and optimum dimensions. In other words, for the range of 
parameters studied (j) is an insensitive parameter within the range studied. 

4. C c /( 1 + e 0 ) is a very sensitive parameter. Up to a certain value of 
the same the cost is not effected much. But, after that with the increase 
of the same cost, increases at a higher rate and beyond a certain value no 
solution can be obtained. 

5. Modulus of elasticity E s has a great effect on optimum design of raft. 
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For very low value of E s it is not possible to make a raft foundation. For 
c — (p -soil when modulus of elasticity value is lower than a limiting value 
and Poisson’s ratio equal to 0.499, solution cannot be obtained and at that 
value of E s , cost of the raft is very high. 

5. The design variable, o ?3 ,i.e., depth of the raft is an insensitive variable 
with respect to all soil parameters taken into account in this study. Even 
if the load is increased final value of d$ remains near about the minimum 
value that is imposed as a constraint. But with the increase of load other 
variable changes and base area of the raft increases. 

6. In case of saturated fine grained soil, for the same column loading 
if the spacing between the two columns decreases then ‘area ratio’ ( A r ) 
increases from the requirement of bearing capacity criteria. But, this not 
so prominent in case of general c — 4> -soil, this is due to the fact that 
imposition of the geometric constraint overrides the design depth required 
from shear and moment consideration. 

3.2 FUTURE SCOPE 

On the basis of the present study one can do the following studies in future. 

1. Optimum cost design of circular and rectangular/square raft treating 
these as plates on elastic foundation. 


2. Optimum design of piled-raft system. 

3. Optimum design of well foundation. 
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